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SUMMARY 

High residual stresses within metal and intermetallic matrix composite systems can 
develop upon cooling from the processing temperature to room temperature due to the 
coefficient of thermal expansion (CTE) mismatch between the fiber and matrix. As a 
result, within certain composite systems, radial, circumferential, and/or longitudinal cracks 
have been observed to form at the fiber-matrix interface region. The compliant layer 
concept (insertion of a compensating interface material between the fiber and matrix) has 
been proposed to reduce or eliminate the residual stress buildup during cooling and thus 
minimize cracking. The present study investigates both elastically and elastic— plastically 
the viability of the proposed compliant layer concept. 

A detailed parametric study was conducted utilizing a unit cell model consisting of 
three concentric cylinders to determine the required character (i.e., thickness and 
mechanical properties) of the compliant layer as well as its applicability. The unknown 
compliant layer mechanical properties were expressed as ratios of the corresponding 
temperature dependent Ti— 24 Al— llNb (a/o) matrix properties. The fiber properties taken 
were those corresponding to SCS-6 (SiC). Results indicate that the compliant layer can 
be used to reduce, if not eliminate, radial and circumferential residual stresses within the 
fiber and matrix and therefore also reduce or eliminate the radial cracking. However, with 
this decrease in in— plane stresses, one obtains an increase in longitudinal stress, thus 
potentially initiating longitudinal cracking. Guidelines are given for the selection of a 
specific compliant material, given a perfectly bonded system. 


INTRODUCTION 

Metal and intermetallic matrix composites are currently being considered for 
advanced aerospace applications due to their attractive high strength— to— density ratio. 
However, due to the inherent coefficient of thermal expansion (CTE) mismatch between 
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fiber and matrix, upon cooling, high residual stresses exist in these composite systems from 
the processing consolidation temperature to room temperature. These residual stresses 
may be large enough in magnitude to generate radial, circumferential (interfacial 
debonding), and/or longitudinal cracks within the matrix. For example, in the case of 
SiC/TijAl+Nb (or alternately Ti-24Al-llNb (a/o)) and SiC/Ti-15-3 systems, 
microscopic radial cracks have been observed to be present at the fiber— matrix interface 
after fabrication and to proliferate after thermal cycling (ref. 1). Figure 1 shows radial 
cracking at the fiber/matrix interface, both after fabrication and after 1000 thermal cycles, 
for the SiC/TigAl+Nb system of interest in this study. In addition to radial cracks, other, 
less frequently observed, crack orientations are illustrated schematically in figure 2(a) 1 . 

A number of potential solutions for reducing the residual stress field have been 
proposed recently. Examples of some potential solutions are high GTE fibers, fiber 
preheating, and the compliant layer concept. The compliant layer concept (the subject of 
the present study) entails the insertion, or addition, of an interface material between the 
fiber and matrix to reduce or eliminate the residual stress field, and therefore the initiation 
of cracks, developed during cooldown. Of particular interest is the reduction in the tensile 
hoop stress a p (fig. 2(b)) within the matrix, which is the cause of the predominantly 
observed radial cracking (figs. 1 and 2(a)). 

Previous investigations have been undertaken to examine the sign and magnitude of 
the resulting residual stress field (refs. 2 to 4) within a composite system as well as the 
effect of including a compliant layer (refs. 5 to 7). Varying degrees of simplification and 
idealization have been adopted, the most common being the assumption of linear elasticity 
with regard to the behavior of the matrix material. In the present study, a variety of 
assumptions with regard to boundary conditions and material behavior are examined 
within the context of a unit cell composed of three concentric cylinders (see fig. 2(b)). 

Here, each cylinder is associated with a different homogeneous isotropic material; that is, 
fiber, compliant layer, and matrix, that may be thermoelastic and/or plastic. The stress 
analysis is a three-dimensional study (assuming perfect bonding between cylinders) and 
considers only the residual stresses due to the initial cooldown cycle during fabrication. 

The objective of this report is to define the required character (i.e., thickness and 
mechanical properties) of a compliant, or compensating, layer which will minimize the local 


! For readers familiar with ASTM E 616—89 standard terminology relating to fracture 
testing, the designation for a radial, circumferential, or longitudinal crack as reported here, 
can alternatively be denoted as (C-R or C-L), (R-C or R-L), or (L-C or L-R), 
depending upon the direction of the crack propagation. 
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tensile residual stresses within the system and thereby reduce the tendency toward 
cracking. The analysis is both analytical and numerical (finite element method) and is 
divided into two parts. The first part contains results associated with a linear elastic and 
thermoelastic stress analysis and parametric study, in which a variety of simplifying 
assumptions are assessed. The second part contains results associated with a 
thermoelastic-plastic finite element stress analysis and parametric study. A section is 
then devoted to identifying candidate compliant layer materials. Finally, the report 
concludes with a discussion of current issues and future work, followed by a summary of 
main conclusions. 


LINEAR ELASTIC STRESS ANALYSIS 

The present investigation begins by revisiting the solution of a linear elastic 
temperature— independent concentric cylinder model both analytically and numerically. 

The motivation behind this starting point is twofold. The first is to obtain an analytical 
solution to study the problem and verify the accuracy (mesh density) of the finite element 
idealization employed; the second is to assess the significance of various assumptions 
imposed by ourselves as well as by previous investigators (refs. 2, 3 and 5); for example, 
boundary conditions, decoupling of in— plane and longitudinal stresses, and 
temperature— independent versus dependent material properties. 

Description of the Analytical Concentric Cylinder Model 

Consider a composite cylinder model, consisting of a single fiber (of radius a) 
embedded in coaxial cylindrical shells of an interface material (outer radius b) and matrix 
material with outer radius c, as shown in figure 2(b). Owing to the obvious cylindrical 
symmetry, we treat the problem in cylindrical polar coordinates r ,0 and z. Assuming the 
application of a uniform temperature change from a processing consolidation temperature 

to room temperature (AT = -1425°F) and planes remaining plane in the z direction, it 
follows from symmetry that the stresses ( cr r , <7 ^, <7 z ) an< ^ strains ( e r i e 0> e z ) are independent of 
the angle ( 6 ) and are functions only of r. Similarly, we assume that each cylinder is 
composed of a homogeneous isotropic material and that a perfect bond between fiber, 
interface, and matrix exists. 

Given these assumptions, the following linear thermoelastic stress— strain relations 
can be written for each cylinder. 
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E 


( 1-2 i/) ( 1 + 1 / ) 


.{ (!— ^)e r + K e ^+ e z ) “ (l+^)aT} 


'* = (l-2,)(l + ,) { + ~ (1+,,)aT> (1) 

a z = K^ r +^) + E(e z - aT) 

where <7 is the Cauchy stress, e is the total strain, v is the Poisson ratio, E is Young’s 
Modulus, a is the coefficient of thermal expansion for the associated material, and T is the 
change in temperature (i.e., T = T — T where T is the current applied temperature and 
T q is the reference temperature). 

Similarly, taking into account symmetry, the required equations of equilibrium 
which the stresses must satisfy are 


£_(* r ) = 0 

dr r 

and equilibrium in the z direction for the entire system is 

r =c 

a x rdr = 0 
z 

= 0 

while compatibility requires that 


/ 


( 2 ) 


( 3 ) 



where the strain-displacement relations are 

e r = i (»W) 

dr 

e - = u ( r ) 

e 0 

r 

and we assume planes remain plane; that is, 

e = constant 
z 


( 4 ) 


( 5 ) 


( 6 ) 
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Note that u(r) is the radial displacement field. 

Upon substitution of equations (1),(2), and (4) to (6), one may obtain expressions 
for the radial displacement and stress components for a generalized plane strain circular 
cylinder (see appendix A). Given these expressions that are valid for a single circular 
cylinder, a variety of laminated concentric cylinder models can be formulated by merely 
applying the appropriate boundary conditions (see table I) and equation (3). 

End and Surface Boundary Conditions 

Figure 3 illustrates schematically the two types of end and surface boundary 
conditions considered in this study; that is, plane strain (e z =0) and generalized plane 
strain (e ^ 0) end constraints and a free surface constraint on the unit cell versus a free 
surface on a homogenized outer cylinder surrounding the unit cell. The rationale behind 
selecting the free surface constraint stems from the thermal nature of the problem, in that 
all the material surrounding the unit cell is undergoing similar thermal strains. The most 
realistic of the two surface constraints is the case when the unit cell is surrounded by an 
effective, homogenized, transversely isotropic, outer cylinder (with a radius sufficiently 
large relative to that of the unit cell), as the influence of neighboring fibers are 
incorporated into the analysis through the rule of mixtures (ref. 8). Results indicate, 
however, that the additional complexity in this idealization is unwarranted since the 
residual stress field developed in the free surface case is sufficiently close to that obtained 
by using the homogenized outer cylinder. This is clearly shown in table II where results 
are given comparing the effect of surface and end boundary conditions as well as a unit cell 
with and without a compliant layer. The geometry (i.e., the radii a, b, and c, representing 

a composite with a 40 percent fiber volume ratio (a /b ) and 10 percent thick (t/a) 
compliant layer, when applicable) and elastic (temperature-independent) mechanical 
material properties for this example are given in table III. 

From an examination of table II, a number of conclusions can be drawn: 

(a) The in-plane stress states (ct. and a 0 ) are basically unaffected by 

modification of the end condition (i.e., e z = 0 or e z = constant) while the 
longitudinal stress (o" z ) is highly affected. In particular, the fiber stress 
becomes compressive under the generalized plane strain condition, while 
under the plane strain condition it is tensile. 
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(b) For a given end condition, modification of the surface constraint results in 
approximately a 13 percent increase in <r r , a 5 to 10 percent decrease in a q, 
and a 20 to 25 percent decrease in a within the compliant and matrix 
material. 

(c) Inclusion of a compliant layer material reduces the tensile hoop stress within 
the matrix (r=b) by approximately 25 percent while a significant stress state 
is produced in the compliant layer. 

Because of the results shown in table II and other results not included here, the generalized 
plane strain end and free surface boundary conditions are assumed when conducting the 
following parametric studies. 


Decoupling of In— plane and Longitudinal Stress Components 

Turning our attention to the results of a previous investigation (ref. 5) wherein the 
in— plane stresses (<r r and < 7 q) were assumed to be decoupled from the longitudinal stress, 
the significance of such an assumption can now be examined, since in the present analysis 
no assumption is made with regard to decoupling. Assuming the geometry and material 
properties to be those taken by Ghosn and Lerch (ref. 5), three concentric unit cell 
calculations have been revisited. The results are tabulated in table IV. 

Upon examining table IV, it is clear that two general observations can be made. 

The first is that the discrepancies in the in-plane stresses (a r and o ^ are less than those in 
the longitudinal component (c^), and secondly, these discrepancies increase as the 
thickness of the compliant layer is increased. Furthermore, if one considers the reduction 
in hoop stress due to the addition of a compliant layer, one sees approximately twice the 
reduction for the cases in which the decoupling assumption is enforced compared with the 
more accurate three-dimensional analysis presented here. Therefore, the impact of this 
decoupling assumption on the resulting stress state is significant and detrimental 
(nonconservative), while its influence on the tractability of the solution is minimal. As an 
aside, a comparison between the effective longitudinal strain that is found by using the rule 
of mixture; that is, 

e z = a P T 


where 


a, = 


V f“f E f + V c a c E c + V m“m E m 


V f E f + V c E c + V m E m 


and that calculated directly from the present analysis (e z ) is included (see table IV) and is 
found to be less than one percent. 
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Description of Finite Element Concentric Cylinder Model 

Four concentric cylinder unit cell models have been idealized utilizing PATRAN 
(ref. 9), each with a different compliant layer thickness; that is, t/a= 0.05, 0.075, 0.1 and 
0.2. Two of these mesh discretizations, those corresponding to t/a = 0.05 and 0.1, are 
illustrated in figures 4 and 5. The element employed is an 8— node isoparametric solid brick 
element with three degrees of freedom per node. A three-dimensional finite element 
analysis, with quarter symmetry 2 was performed using the finite element program MARC 
(ref. 10). The end and surface boundary conditions assumed were those of generalized 
plane strain and free, respectively, as discussed previously. Due to the symmetry of the 
problem, the natural coordinate system of choice is that of polar coordinates; therefore, 
transformation of the resulting stress and strain fields from the assumed (prescribed) global 
cartesian coordinate system of MARC is required. The resulting principle stress and strain 
components are then processed further by calculating the arithmetic mean between two 
integration points within a given element (see fig. 6). Note that an eight— point Gauss 
quadrature numerical integration scheme was selected that yielded four data points of 

interest per element, two within a given z plane. 

A mesh density study was conducted to verify the accuracy of the results. The final 
mesh discretizations employed are those shown in figures 4 and 5. An illustration of the 
discretization accuracy (less than one percent in error when compared with the exact 
solution) is shown in figures 7, 8, and 9 in which the radial, circumferential, and 
longitudinal stresses are displayed as a function of radius for the case described in table III 
and idealized in figure 5. 

Temperature-Independent Versus Dependent Material Parameters 

Previously, only temperature— independent material properties have been considered 
for each constituent. Experiments indicate, however, that the material properties are in 
fact temperature— dependent, that is to say, the stiffness and CTE of the matrix and the 
CTE of the fiber vary with temperature. Table V summarizes the temperature-dependent 
behavior of the TijAl+Nb matrix (elastic and plastic properties) and that of the SiC fiber 
(elastic only). These material data were obtained through private communications with 
Brindley (ref. 11) and Textron Inc (ref. 12). In order to assess the importance (with regard 
to modifying the resulting unit cell stress distribution) of incorporating temperature- 


2 The current problem actually possesses full rotational symmetry (i.e., is 
independent of the angle 0), thereby allowing the alternative use of axisymmetric elements. 
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dependent material behavior in the analysis, the problem with parameters given in table III 
was reanalyzed assuming both temperature— independent and dependent material 
properties. Note, as before, that the compliant layer properties are taken relative to those 
of the matrix for all temperatures. 

Table VI contains the results of the analysis. As indicated, the difference between 
the temperature-dependent and independent results is approximately 11 percent. 

However, this difference would increase as the temperature dependence of the material 
properties increased. This can be illustrated theoretically by writing the incremental 
stress-strain relations (see appendix B). On looking ahead to conducting an elastic— plastic 
analysis, it is clear from Table V that the plastic material properties are highly 
temperature dependent. Therefore, material temperature dependence has been included in 
all remaining analysis results. 

Compliant Layer Thermo— Elastic Parametric Study 

Thus far only an assessment and validation of fundamental solution assumptions 
have been examined. The questions remain as to the impact of the so-called compliant 
layer on reducing the overall residual stress field within a composite system and as to its 
required character (i.e., thickness and mechanical properties). 

Interpretation of Problem 

Let us begin by considering, generically, the interaction of the fiber and matrix in 
the r — 6 plane (see fig. 10(a)) and let us assume, as is the case for the present system of 
interest (i.e., SiC/TigAl+Nb) that the CTE of the matrix is greater than that of the fiber. 

In this figure, a free body diagram is shown in which the initial fiber/matrix system 
present at the consolidation temperature T Q is disassembled. The location of each 
constituent, at the reference temperature, is denoted by the solid lines. Upon cooling, both 
the fiber and matrix will contract (assuming their CTE’s to be positive), with the matrix’s 
contraction being greater than that of the fiber by an amount 6. 

Superposition of fiber and matrix (by enforcing compatibility and thus eliminating 
this S overlap) results in the development of stresses within the system. These stresses are 
such that the radial and circumferential fiber stresses are compressive, and the radial and 
circumferential matrix stresses are compressive and tensile, respectively. Consequently, 
the potential for the formation of radial cracks (fig. 2) exists due to the tensile hoop stress 
within the matrix. Whether or not cracking will occur clearly depends upon the 
magnitude of 6 (CTE mismatch) and the ductility and strength of the matrix, or more 
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precisely, the fiber/matrix interface, because the maximum circumferential stress occurs at 
the interface and diminishes with the square of the radius (cf. equation A3 in appendix A). 

Ideally, this problem could be avoided by the insertion of a layer of air 6 thick 
between the fiber and matrix, thereby compensating for the mismatch in thermal strains 
attained during cooldown (see fig. 10(b)). In other words, the development of thermal 
stresses within the system can be mitigated by allowing unrestrained thermal strains to 
occur. Obviously, the value and manufacturing ability of such a composite system is 
unrealistic. However, this ideal situation does suggest that insertion of a material with the 
"proper character" between the fiber and matrix should minimize, if not eliminate, the 
buildup of thermal stresses within the plane (see fig. 10(c)). 

As previously stated, the question at hand is the determination of this "proper 
character." Initially, if we restrict our attention to thermoelastic behavior, this character 
can be described mechanically by three parameters: the coefficient of thermal expansion, 
the thickness of the layer, and the stiffness. Utilization of figure 10 can provide insight 
with regard to the preferred magnitude of the unknown compliant layer’s coefficient of 

c c 

thermal expansion (a ), assuming a given thickness tfplus and stiffness E . 

For example, if a c is equal to a (the CTE of the fiber), a similar if not identical 
resistance to the contractive thermal strains within the matrix exists as in the system with 

c 

no compliant layer. Alternatively, if the unknown layer were to possess an a such that it 
was between that of the fiber and matrix, some reduction in thermal strain mismatch 
(displacement 6) could be achieved. However, if the CTE of the unknown layer were 
greater than that of the matrix, the contraction of this layer would be greater than that of 
the matrix and thereby allow the matrix to contract freely. Restated, the unknown 
material would in essence "get out of the way" of the surrounding matrix material and thus 
resemble the ideal case in which no thermal stress is generated. Clearly , due to 

compatibility constraints, if a c is sufficiently large and perfect bonding exists between the 
compliant layer and matrix, the matrix may indeed be pulled, or forced, to contract beyond 

the usual a m AT amount, thus inducing a tensile radial and compressive hoop stress state 

within the matrix. Such an argument would suggest that there exists an a (greater than 

a m ) and thickness combination which would eliminate the in— plane stress state within the 
matrix while increasing the stress in the fiber. Intuition also tells us that the stress state 

c c 

within the compliant layer will increase as a increases, given a fixed thickness and E , 

since in essence the load on the compliant layer is being increased. 
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Parametric Study 

A parametric study was conducted involving the variation of the CTE, thickness, 
and modulus to substantiate the above observations relative to thermoelastic behavior. 
Selected results are presented in figures 11 to 14. Figures 11 and 12 illustrate the 
variation of circumferential and longitudinal residual stress (i.e., the stress attained upon 
reaching room temperature) as a function of the ratio of the CTE of the compliant layer 

and matrix (i.e., a c /a m ), respectively. This CTE ratio can be extended to incorporate 
other fiber— reinforced systems by utilizing an alternative normalized CTE ratio; that is, 

(a c -a)/(a m -a), wherein the driving a-mismatch between that of the fiber and matrix is 
incorporated. 

The validity of the stresses shown correspond to the inner radius r^a and r^b for the 
matrix, respectively, while those associated with the fiber are valid for all locations within 
the fiber. Also included in the figures are the results for two normalized compliant layer 
thicknesses, t/a = 0.1 and 0.05, where t is the thickness of the compliant layer and a is the 
radius of the fiber. Note that the stresses reported in figures 11 and 12 are associated with 

a compliant layer whose modulus (E c =0.5E m ) and Poisson ratio (i/ c =v m ) are held fixed. 

A comparison of these stress curves within any given material is not strictly valid 
since the location of the corresponding integration points, the inner radius, changes with a 
variation in compliant layer thickness. This variation is quite small, however, and for the 
purposes of this study such a comparison is justifiable. 

An examination of figure 11 shows that increasing the normalized CTE ratio for a 
given thickness, t/a, decreases the circumferential stress in the matrix while increasing it in 
the compliant layer and fiber 3 . Furthermore, figure 11 illustrates that, for a given CTE 

ratio (e.g., a c /a m = 2.0), increasing the thickness, t/a, of the compliant layer results in a 
decrease in the magnitude of the circumferential stress within the matrix, compliant layer, 
and fiber. Similar (but less significant) trends with regard to the change in longitudinal 
stress are observed in figure 12. Figures 11 and 12 both illustrate a rapid increase in stress 

within the compliant layer as its CTE, a c , is increased. Experience suggests that this 

increase could be compensated for by modifying the elastic stiffness, E , of the compliant 
layer. 


3 Note that a CTE ratio of 1.0 corresponds to the case when the CTE of the 
compliant layer and that of the matrix is the same, while a ratio equal to 0.392 corresponds 
to the case when the compliant layer’s CTE is equal to that of the fiber. 
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Figures 13 and 14 illustrate the effect of varying the stiffness on the circumferential 
stress (crfj) distribution when a constant thickness, t/a=0.1, is specified. The effect on the 
radial distribution of Cq is summarized by plotting the hoop stress induced at both the 
inner radius, IR, and outer radius, OR, of the compliant layer cylinder (fig. 14). Two 
dominant trends are exhibited in figures 13 and 14. The first is that the magnitude of the 
circumferential stress for a given CTE ratio is significantly decreased within the compliant 
layer as well as within the fiber when the stiffness of the compliant layer is decreased, while 

the matrix remains basically unaffected for a c / a m < 4. The second trend deals with the 

increasing influence that stiffness, E c , has on the resulting stress distribution (i.e., the 
difference between IR and OR stress magnitudes increases as stiffness increases) when the 
CTE of the compliant layer is increased. 

Experimental data indicate that the matrix’s yield stress is approximately 50 ksi at 
room temperature. It is therefore clear from the high stress values computed (those values 
greater than 50 ksi) in the compliant layer and matrix material that a purely thermoelastic 
analysis is insufficient to provide realistic results. Therefore, the next logical step is to 
incorporate yielding, or plasticity, of both the unknown compliant layer and the matrix 
material into the parametric stress analysis. 

Summary of Thermoelastic Trends 

A simple thought experiment was described to help explain and motivate the use of 
a larger CTE in the compliant layer than that in the matrix. Calculations showed that: 

(1) Increasing the normalized CTE ratio of the compliant layer decreases the 
circumferential and longitudinal stress in the matrix while increasing it in 
the compliant layer and fiber. 

(2) Increasing the thickness (t/a) of the compliant layer results in a decrease in 
the magnitude of circumferential and longitudinal stress within the matrix, 
compliant layer, and fiber. The longitudinal fiber stress is, however, 

increased beyond an a c / a m ratio of two. 

(3) Decreasing the stiffness of the compliant layer has a significant impact on 
reducing the circumferential stress in the compliant layer and fiber, while the 
matrix stress is relatively unaffected. 

(4) The influence of the stiffness increases as the CTE of the compliant layer is 
increased. 

The reader is cautioned not to draw any definitive conclusions at this point since 
some important departures from the above trends are observed when plasticity is 
incorporated in the analysis, as illustrated in the next section. 
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THERMOELASTIC PLASTIC STRESS ANALYSIS 

The incorporation of yielding, or plasticity, is accomplished by assuming a uniaxial 
bilinear stress— strain response of the matrix, as shown in figure 15 and table V. Note that 

the yield point denoted by is defined by using a 0.02% offset criterion, and H m defines 
the linear hardening slope. The temperature-dependent material response determined 
experimentally by Brindley (refs. 11 and 13) for the Ti^Al+Nb matrix clearly supports this 
assumption. Again the unknown compliant layer material properties are taken relative to 

c m c in 

that of the known matrix material throughout the temperature range E /E o / q 1 , and 

H c /H m . The multiaxial yield criterion utilized was that of von Mises. All stress analyses 
were performed numerically with the nonlinear finite element program MARC 4 

The inclusion of plasticity effects within the compliant layer and matrix were found 
to have a significant impact on the magnitude and distribution of stress within the 
composite system. This impact is best illustrated by example. For comparison purposes 
both the elastic and elastic— plastic calculations are shown in table VII, in which the 
distribution of radial, circumferential, and longitudinal stress versus radial location are 

tabulated, for the compliant layer idealization depicted in figure 5 with a c / a m = 2.0, 

E c /E m =0.5, c7y/a^ l =0.5, and H c /H m =0.5. An examination of table VII reveals that 

inclusion of plasticity in the calculations results in approximately a 20 to 70 percent 

reduction in stress magnitudes when compared with the purely elastic calculations. As one 

might expect, the largest reductions are confined to areas of highest stress, or largest 

plastic strains, as in the compliant layer. Another point of interest is the redistribution of 

the longitudinal stress a . In the elastic case, this stress is constant throughout a given 

material, while in the plastic case the stress drops at the inner radius due to the drop in 

radial and circumferential stress (see equation (1)) and then increases as one moves toward 

the outer radius. Note that in the matrix cylinder, a approaches that value associated 

with a purely elastic solution. This is due to the fact that the outer radial region remains 

elastic until the final increment in temperature drop. Thus if AT were larger or the matrix 

softer, then a at r=c would be lower, as illustrated later, 
z 


4 A simplified analytical solution can be obtained if one assumes a Tresca yield criterion and 

assumes that a < a < a* However, except for a small region near the inner radius of the 
r z c 1 

compliant layer (table VII), this latter assumption is invalid for the majority of cases under 
investigation. Therefore, using the Tresca criterion does not provide the analyst with any 
greater analytical simplicity than use of the von Mises criterion. 
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Compliant Layer Thermoelastic Plastic Parametric Study 

In light of including material nonlinearity, let us again examine the role of a 
compliant layer in reducing the overall residual stress field within a composite system, in 
particular, the tensile hoop stress which is believed to be linked to the radial cracking of 
the fiber-matrix interface. The problem definition discussed earlier and shown in figure 10 
still applies and suggests that the CTE and thickness of the compliant layer will play the 
dominant role with regard to reducing the tensile hoop stress within the matrix. 

This was investigated by varying the CTE, a c , and thickness, t/a, of the compliant 

c c c 

layer while holding the remaining available parameters E , o , and H fixed. These 

V 

parameters were taken to be 0.5 relative to those of the matrix throughout the temperature 

range E c /E m =0.5, a c /a m = 0.5, and H c /H m = 0.5. 
y y 

Figures 16, 17, and 18 illustrate the variation of circumferential, longitudinal, and 
radial stress components for each constituent as a function of the normalized CTE ratio 
c in 

a / a , respectively. The stress values shown correspond to the inner radius for the 
compliant and matrix cylinder, while those associated with the fiber are valid for all radial 
locations within the fiber. Furthermore, results associated with four normalized compliant 
layer thicknesses (t/a = 0.05, 0.075, 0.1, and 0.2) are included in figures 16 to 18. 

One can see from figure 16 that, for a given normalized thickness t/a, increasing 
the normalized CTE decreases the circumferential stress in the matrix and fiber while 
increasing the stress in the compliant layer. Notice that beyond a given CTE ratio (e.g., 

a c /a m — 2.0 for t/a = 0.1) the compliant layer hoop stress becomes greater than that of 
the matrix. The actual magnitude of this compliant layer stress component is immaterial 
at this stage, since with proper material selection this increase in stress can be controlled. 
This will become clear in the following section. The second major trend observed is that 
given a CTE ratio greater than 0.392, increasing the normalized thickness decreases the 
matrix and fiber hoop stress while increasing the compliant layer stress. As alluded to 

earlier, there appears to exist an a and thickness combination that can eliminate the 
tensile hoop stress in the matrix. Therefore, the net result of including a compliant layer is 
to lower, if not eliminate or change, the sign of the tensile hoop stress component of the 
matrix and fiber. Clearly this reduction is at the cost of increasing the compliant layer 
stress. 
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A comparison of figure 16 with its purely elastic counterpart (fig. 11) shows a 
number of differences. The first is the nonlinearity of the stress versus CTE curves 
corresponding to the elastic/plastic case when compared with those of the elastic. The 
second is that the fiber stress in the elastic cases increases with respect to increasing CTE, 
while the fiber stress corresponding to the plastic cases decreases. This reversal in trends 

can be explained by considering a fiber subjected to a pressure loading equal to P c = 
c c • c 

E o AT. Thus, if a is increased while holding the stiffness and change in temperature 
constant, one would expect the stress in the fiber to increase as in figure 11. However, 

when plasticity is included, this pressure P c is limited by the yielding of the material and 
in fact, due to the coupled kinematic effects, this pressure is actually reduced along with 
the stress within the fiber. The third difference is a subtle shift in the thickness 
dependence in the compliant layer stress. The fourth difference, as noted earlier, is an 
overall lowering of the stress values with the inclusion of plasticity. The general trend, 
however, remains intact with regard to reducing the matrix stress by increasing either the 
CTE or thickness of the compliant layer. 

In examining the longitudinal stress versus normalized CTE ratio (fig. 17), one 
observes within all three constituent materials an increase in longitudinal stress with an 
increase in CTE ratio. The longitudinal stress in the matrix and fiber does, however, 
reverse this trend beyond a certain critical CTE ratio, determined by the change in sign of 
the sum of the circumferential and radial stress. Similarly, an increase in longitudinal 
stress is observed with an increase in normalized thickness (for a specified CTE ratio) in 
the matrix and compliant constituent materials. Note that increasing the thickness lowers 
the critical CTE at which the longitudinal stress in the matrix peaks out. 

In comparing the trends indicated in figures 12 and 17, a number of differences are 
observed that are similar to those noted earlier for the circumferential stress. These 
include the nonlinear appearance of the curves in the plasticity case, the reversal in fiber 
stress with respect to increasing the CTE, a reversal in thickness dependence, and an 
overall reduction in stress levels with the inclusion of plasticity. The most important 
discrepancy between an elastic (fig. 12) and plastic (fig. 17) analysis resides in the trends in 
the matrix. In the elastic case, increasing the CTE ratio or normalized thickness decreases 
the longitudinal stress, whereas when plasticity is included, the stress is increased. This 
increases the potential for initiating cracking in the longitudinal direction. Provided 
loading is in the fiber direction, longitudinal cracking is considered to be more detrimental 
then radial cracking, thereby requiring that upper limits on both thickness and CTE of the 
compliant layer be imposed. These limits will be discussed later. 
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Figure 18 illustrates the effect on the radial stress distribution of varying the CTE 
and thickness of the compliant layer. Clearly, as the CTE and/or thickness are increased, 
the radial stress throughout the system is reduced. A comparison of trends resulting from 
an elastic (fig. 11) analysis with those from the plastic analysis reveals a nonlinear 
appearance in the curves for the plasticity case and a reduction, instead of an increase, in 
radial stress with respect to increasing CTE. This reduction in radial stress with increasing 
CTE suggests another potential problem regarding the use of compliant layers; that is, 
bond strength. 

Throughout this study an assumption of perfect bonding has been made. However, 
in actuality, a wide range of bonding strengths can exist. For instance, if the bonding 
strength is purely mechanical (e.g., frictional), reducing the radial stress is proportional to 
reducing the bond strength, and therefore the load carrying capability, of the composite 
system. As a result, an upper limit on the allowable CTE and/or thickness of the 
compliant layer may be necessary if the bond strength at either interface is questionable. 
For example, if the bond strength is only frictional, then this limit would be reached when 
a the radial stress, becomes zero or tensile. 

Effect of Varying Compliant Layer Material Parameters 

In the previous section only variations in compliant layer CTE and thickness have 
been addressed, with all other parameters held fixed. As observed previously, increasing 
the CTE or thickness resulted in an increase in compliant layer stress in the longitudinal 
and circumferential directions. Therefore, the question at hand is what impact do the 
material properties of the compliant layer have on the resulting residual stress distribution. 
To answer this question, seventeen different sets of material parameters were investigated. 
Cases in which the stiffness, yield point, and hardening slope are all normalized with 
respect to the corresponding matrix properties are described in table VIII. Results 
indicating the radial, circumferential, and longitudinal stress at specific radial locations are 
shown for all seventeen cases in tables IX to XI. These results are associated with a 

specific normalized thickness, t/a = 0.1, and CTE ratio, a c /a m = 2.0. 

It is apparent from these tables that variations in mechanical properties have little 
if any impact on the fiber (roughly 20 percent) and matrix (under 5 percent) stress state, 
while the compliant layer stress state is greatly (as much as 90 percent) affected. 

Therefore only the compliant layer stress state as a function of material property will be 
studied. The variation in circumferential and longitudinal compliant layer stress versus 
normalized hardening slope, yield point, and stiffness is illustrated in figures 19 to 21. 
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In figure 19, the effect of varying the hardening slope H c /H m (between 0 and 1) on 

the circumferential and longitudinal stress is shown, given a relative yield point a y/ 

= 0.5 and stiffness E c /E m = 0.666. As expected, the perfectly plastic case (i.e., H c /H m = 
0 ) gives the lowest stress state, while an increase in the hardening slope increases the 
stress state. Note also that increasing the hardening slope causes an increasing stress 
distribution; that is, the difference between the stress at the inner and outer radius 
increases. 

Figure 20 illustrates the effect of varying the yield point on the 

c m 

circumferential and longitudinal stress distribution for a given stiffness E /E = 0.666 and 

hardening slope H c /H m = 0.5. Again as one might suspect, the stress state is reduced as 
the yield point is reduced. The deviation between inner and outer radial stress values 
decreases with increasing yield for the circumferential component, while it increases with 
increasing yield for the longitudinal component. 

Figure 21 illustrates the effect of varying the stiffness E c /E m on the compliant layer 
stress state, given a relative yield and hardening slope of 0.5. Once again it is observed 
that increasing the stiffness ratio increases (at a decreasing rate) the circumferential and 
longitudinal stress components. Note that only the value of stress, not the distribution, is 
affected by modification of the stiffness ratio, and that the percent increase or decrease is 
much less per increment in stiffness ratio than that of the yield or hardening slope. 

To gain insight into the relative importance of each available material parameter 
with regard to reducing the stress state within the compliant layer, the circumferential and 
longitudinal stress components are put in ascending order, from minimum to maximum 
(table XII). Note that the material parameters corresponding to case number 16 produce 
the minimum stress state, both circumferentially and longitudinally, with respect to inner 
and outer radial location, while case number 2 produces the maximum stress values. Table 
XII also contains the ratio of material properties associated with each case number, similar 
to table VIII but listed in ascending order of the resulting stress values. Column 1 
represents a measure of the interaction between the hardening slope and yield point of the 
material, and is associated with the sum of the two ratios. Clearly, the smaller this sum 
the lower the stress state within the compliant layer. Columns 2, 3, and 4 are associated 
with the relative yield point, hardening slope, and stiffness ratios, respectively, and, as in 
column 1, the lower the ratio (provided the other two ratios remain fixed) the lower the 
stress state within the compliant layer. 
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With regard to impact on the resulting residual stress state, the level of importance 
increases as one moves from the right column to the left. This level of importance is 
illustrated by considering cases in which columns 1 and 4 are the same while the yield and 
hardening ratios of columns 2 and 3 vary. Under these constraints, cases exist which 
indicate as one might expect that more importance should be placed upon having a lower 
yield point than upon having a lower hardening slope (see cases 1, 10, 14, and 15 in table 
XII). The inherent coupling between yield point and hardening slope is illustrated by cases 
16, 4, and 11. This prompted the inclusion of a still more important parameter, namely, 

the sum of the yield point and hardening slope. 

When examining cases 10 and 12, one might observe that an exception to the above 
level of importance is indicated in table XII. Here it is not the lower yield point that 
dictates the lower stress state but rather the hardening slope. As noted in table XII, 
however, this is true only for the circumferential and longitudinal stress components 
associated with the inner radius; those associated with the outer radius follow the above 
guidelines. 

Summary of Trends 

In review, calculations have shown that the CTE and thickness dominate when 
attempting to reduce in-plane fiber and matrix stress. More specifically, increasing the 
CTE and thickness of the compliant layer decreases the in-plane stress within the 

fiber and matrix while increasing the out— of— plane (& z ) stress component. The compliant 

layer’s stiffness E c , yield point 0 ^, and hardening slope H c dominate when attempting to 
reduce the stress state within the compliant layer, and yet have little if any effect on the 
fiber or matrix stress state. Furthermore, it has been observed that, when selecting a 
compliant layer material, the properties should be chosen such that 

(1) a > a 

(2) t /a should be as large as other considerations allow 

(3) H c /H m + a c fa m should be small 

y y 

(4) Yield point a £ relative to matrix should be low 

(5) Hardening slope H c relative to matrix should be low 

(6) Elastic stiffness E c relative to matrix should be low 

The above list is in order of importance, with respect to impact, for obtaining a minimum 
overall residual stress state. Differences between performing elastic and elastic— plastic 
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analyses have also been noted. These include 

(1) Nonlinear appearance of the curves in the case of plasticity, 

(2) Reversal in fiber stress and longitudinal matrix stress with respect to 
increasing the CTE, and 

(3) Reversal in thickness dependence and an overall reduction in stress levels 
with the inclusion of plasticity. 

MATERIAL SELECTION 

Given the above material parameter guidelines, candidate compliant layer and/or 
diffusion barrier materials can be put forth. Prior to selecting these candidate materials, 

more quantitative bounds on the key parameters a c /a ra > 1.2 and 0.1 < t/a < 0.2 were 
obtained by calibrating the analysis with a compliant layer to one without a compliant 
layer (i.e., SiC/TigAl+Nb). 

Results of the fiber/matrix analysis indicated that the circumferential stress at the 
fiber/matrix interface was 44.5 ksi, while the longitudinal stress at the inner and outer 
radius of the matrix cylinder was 39 and 62 ksi, respectively. Experimental results indicate 
that radial cracks formed at the fiber-matrix interface in the SiC/Tig A1 + Nb system, 
while few if any longitudinal cracks appeared. As a result, the present analysis is 
calibrated to the experimental observations by assuming that the fracture strength at the 
inner diameter of the matrix is equal to 44.5 ksi, while that at the outer diameter is greater 
than 62 ksi. In other words, a gradient in fracture strength is applied to the simplified 
analysis to account for the unknown interface properties as well as possible material defects 
in the actual system. Thus a circumferential stress greater than or equal to 44.5 would 
result in radial cracking, while a longitudinal stress less than 62 ksi would not result in 
axial cracking. The lower limit on the CTE ratio was determined by imposing this limit 
of c q < 44.5 on the curves of figure 16. Similarly, imposing the limit that a z < 62 on the 
curves of figure 17 gives the upper limit on the thickness. Note that the stress values given 
in figure 17 correspond to only those associated with the inner radial locations, although as 
shown in figure 22 the stress associated with the outer radial location corresponds to those 
of the inner just prior to the peaking of the longitudinal stress. 

These bounds should not be considered absolutes, since the interface material in the 
actual SiC/Tig Al+Nb composite is unknown (thus critically affecting the stress values) 
and the given stress values are at room temperature (cracking could have initiated at a 
higher temperature and therefore lower stress). Furthermore, a lower upper bound on the 
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compliant layer thickness can be achieved based on overall composite density 
considerations. 

Considering the above limits and the fact that processing temperatures are typically 

above 1500°F (thus requiring a melting point in excess of 1500°F), pure element and alloy 
candidate compliant/compensating layer materials are put forth in tables XIII and XIV, 
respectively. These candidates are limited still further by practical considerations such as 
handleability in air, radioactive elements, carcinogenic properties, and potentially too high 
a yield to maintain adequate stress levels in the compliant layer itself. Stars in the above 
tables indicate the affected candidate materials. With respect to pure elements, the 
candidates of choice appear to be silver, copper, and gold, with copper being the best of the 
three, in the authors’ opinion, from a mechanics and economics standpoint. 

Chemical compatibility with the fiber and matrix is probably the most restrictive 
criterion for a candidate compliant layer material. Thermodynamic studies (ref. 14) 
conducted on various candidate materials have revealed that copper is indeed the system of 
choice provided a diffusion barrier between the copper and Ti^Al+Nb matrix is used. The 
diffusion barrier selected for the present study is niobium. This system will be fabricated 
by Textron and tested by Lewis personnel in order to experimentally verify the feasibility 
of the compliant layer concept. 

DISCUSSION OF ISSUES AND FUTURE WORK 

This report has taken us one step closer to understanding the advantages and 
disadvantages of employing a compliant /compensating layer. A number of issues still 
remain to be addressed. First is the effect of fiber interaction, spacing, and volume 
fraction. Second is the effect of warping (i.e., the removal of the assumption that planes 
remain plane) and fiber breakage. These must be examined since such a study would 
introduce longitudinal shears into the stress analysis and allow us to address the issue of 
circumferential debonding and load transfer (increase in critical length) caused by the 
inclusion of a compliant layer. Third is the effect of multiple compliant layers, since all 
conclusions stated in this report are valid only for a single layer. Fourth, the stress 
analysis must account for material time dependency during the fabrication cycle in order to 
accurately determine the true residual stress field and draw conclusions about optimum 
cooldown histories. This would require the inclusion of viscoplasticity in the stress 
analysis. Fifth, we must address the question of the fabricability and effectiveness of 
utilizing a thermally anisotropic interface layer whereby the radial and longitudinal CTE 
are different. The assumption here is that if the longitudinal CTE of the layer were 
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between that of the fiber and matrix, and the radial CTE were greater than the matrix, 
then the hoop stress within the matrix could be reduced while not adversely affecting the 
longitudinal stress. Lastly, prior to passing final judgment on the adequacy or inadequacy 
of utilizing a compliant layer, thermal cycling with and without a mechanical load needs to 
be addressed. This will be attempted by using the concept of shakedown. All of these 
areas are presently being investigated. 

SUMMARY AND CONCLUSIONS 

The present study has investigated both analytically and numerically the use of a 
compliant layer (or better named, compensating layer) concept to reduce, if not eliminate, 
residual stress buildup, and therefore cracking, during fabrication cooling. Significant 
effort has been devoted to understanding and describing the mechanics of the problem, 
specifically, the importance of the initial stress analysis assumptions. A detailed 
parametric study was performed using a finite element concentric cylinder model with 
generalized plane strain end conditions and free boundary conditions. The fiber, SiC, was 
assumed to be isotropic and linear elastic, while the unknown compliant layer and matrix 
(Ti^Al+Nb) cylinders were assumed to be isotropic and bilinear elastic— plastic. The 
unknown compliant layer properties were taken to be relative to the corresponding known 
matrix properties over the full temperature range. Finally, perfect bonding between 
cylinders was also assumed and only the initial cooldown cycle was analyzed. Within the 
confines of these assumptions, the following general conclusions have been reached: 

1. Plasticity was shown to have significant influence on the actual magnitude 
and distribution of stress (i.e., the resulting trends) and therefore must be 
included in the analysis. 

2. A compliant/compensating layer can be used to reduce, if not eliminate, 

in— plane residual stresses within the fiber and matrix and therefore 

radial cracking as well. 

3. With this decrease in in— plane stresses comes an increase in longitudinal 
stress, thus potentially initiating longitudinal cracking. 

4. Similarly, if bonding was strictly mechanical, especially frictional, the 
reduction in radial stress with the addition of a compliant layer could be 
extremely detrimental to the overall composite performance. 
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5. The two most important interface layer parameters for reducing the matrix 
in-plane stresses are the CTE and thickness of the layer; thus the reason for 
the more appropriate name — compensating layer. These two parameters 

should be chosen such that a c > a m and the thickness, t/a, is as large as 
other considerations will allow. 

6. The mechanical properties (yield point, hardening slope, and stiffness) play 
the dominant role in reducing the stress state within the interface layer and 
should be chosen such that relative to the matrix they are as low as possible. 

7. When selecting an interface material, the order of importance, with respect 
to impact, for obtaining a minimum overall residual stress state, is, 

a c > a™, t/a, H c /H m + ape™ , a*, H c , and E c . Be aware, however, that 
this requirement may not provide maximum life under cyclic conditions. 

8. A candidate system (SiC/Cu/Nb/Ti3Al+Nb) has been selected to be 
manufactured and tested experimentally to verify the applicability of 
employing a compliant/compensating layer technology. 
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APPENDIX A 


By substituting equations (1), (2) and (4) to (6), one can obtain the general 
equations describing the radial displacement, u, and principal stress components, cr r , a ^ 

a , for a circular cylinder with generalized plane strain end conditions. These expressions 
z 

are 


ufr) = ( 1 + u ) a f Trdr + C,r + C 2 

TT^JT\ 1 7 


(Al) 


a r ( r ) 


-E a f Trdr + E [( C l + 

( l-i/ ) r 2* / r. TTT7)[ (1-2 v) 


( C 1 + ye z) - CJ 1 ) 


(A2) 


a 




E g f Trdr- EaT 
( l-*') r 2 rj ( l-v) 


E 


( C 1 + ye « ) + C„(J_) 


( l + v ) L (1—2 v ) 


(A3) 


a — 


— E aT 


+ 


E 


( 1 + v) ( 1-2 v ) ( l + v)[ 


2vC 1 + (1— i/)e 2 


(A4) 


With the above expressions valid for a single circular cylinder, a variety of different 
laminated concentric cylinder models can be formulated by merely applying the 
appropriate boundary conditions to determine the constants of integration (C^ and Cg). 

Consider, for example, the three concentric cylinder case shown in figure 3(a) in 
which the generalized plane strain end and free surface boundary conditions are assumed. 
In this case, six unknown constants of integration (Cj and C 2 corresponding to each 

constituent, i.e., fiber, compliant, and matrix cylinder) and the axial strain e z must be 

determined before the radial displacement and three principle stress distributions will be 
known throughout the unit cell. Thus the seven boundary and compatibility conditions 
listed in table I must be satisfied. 


Applying conditions (a) through (g) of table I to equations (Al), (A2), and (A4) 
results in the following seven constraint equations: 

C* = 0 (a) 


(C?-Cf)a 2 + C$ = ( 1 + I/ f) or f^rdr 

11 2 TT^fT { J 

1 n 


0 >) 
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Writing the nontrivial equations (equations (b) to (g)) in matrix form results in: 

[M] {X} = {R} . 

which can be solved by a simple inversion process; that is, 

{X} = [M]-‘{R} 

where the vector {X} T = {c[ e z ) contains the constants ot integration for 

equations (Al) to (A4). The solution to this system, as well as the displacement and stress 
field within each cylinder, has been accomplished by using the software package 
MATHCAD (ref. 15). The MATHCAD document containing the above details is given 
below. Similar boundary conditions can be applied to obtain displacement and stress 
distributions for the other boundary conditions discussed in this report. 


MATHCAD DOCUMENT 


ANALYTICAL SOLUTION FOR CONCENTRIC CYLINDER MODEL - THREE MATERIALS 
INNER: FIBER; MIDDLE: COMPLIANT; OUTER: MATRIX 

Generalized plane strain assumption 

Free surface constraints 

DEFINE MATERIAL PROPERTIES AT ROOM TEMPERATURE; REFERENCE TEMPERATURE 

IS AT 1500 F 


FIBER MATRIX COMPLIANT LAYER 


EF 

:= 58000. 

EM := 16000. 

RE 

:= 0.5 


-6 

-6 



AF 

:= 1.96- 10 

AM := 5.0- 10 

RA 

:= 2.0 

VF 

:= .25 

VM := .26 

VC 

:= 0.26 


EC := EM- RE 
AC := RA-AM 


DEFINE GEOMETRY AND TEMPERATURE CHANGE 


fiber OD compliant OD matrix OD 

a := 0.632 b := 0.6952 c := 1.0 


change in temp absolute DT 

DT := -1425. ABDT := 1425 
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Fiber 


SOLUTION OF CONSTANTS OF INTEGRATION 

Matrix Compliant layer 


EF 


A := AF- 


1 - VF 


EM 


C := AM- 


1 - VM 


EC 


CC := AC 


1 - VC 


EF 


B : = 


EM 


D : = 


1 + VM 


(1 + VF) • (1 - 2- VF) 

DM 


EM 


EC 


CD : = 


1 + VC 


EC 


(1 + VM) • (1 - 2 VM) 


CDM : = 


(1 + VC) • (1 - 


Longitudinal strain (left side) 

2 2 2 2 2 
a b - a c - b 

EZ6 := B (1 - VF) + CDM (1 - VC) + EM- (1 - VM) 

2 2 2 


Longitudinal strain (right side) 

2 [2 2 ] [22 

REZ := A a + CC- |b - a J + C- [c - b 


APPLY BOUNDARY CONDITIONS TO FORM MATRIX AND OBTAIN 
APPROPRIATE COEFFICIENTS 


M2 6 


'-B VC 

VF + 

CD 1 - 2 - VC 


M3 6 


-CDM VM 

VC + 

D 1 - 2 - VM 


2 - VC) 
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2 


M : = 


-a 


-B 2 

a 

_-i 

— a 
CD 

0 

1 - 2 - VC 
2 
-b 

-1 


-CDM 2 

- w 

CD 

0 

JJ 

D 

D 

0 

0 

0 


1 - 2VM 
2 
c 

1 - 2VM 


-1 


B VF- a CDM- VC- lb - a 


2 2] r 2 2 1 

b - a J 0 DM-VM- |c - b J 


M26 a 


-1 M3 6- b 


VM 

1 - 2- VM 
EZ6 


define the right-hand side 


DT2 : = 


DT 

2.0 


R : = 


i + vf ; 

AF DT2 • a 

1 - VF 

-A 2 

— DT2 ■ a 
CD 

1 + VC f 2 21 

AC- DT2- b - a 


1 - VC 


f ■ .1 

-CC [2 2] 

DT2 • |b - a J 

[= 2 - A 


D 

1 + VM r 2 2 

AM DT2 I c - b 


1 - VM 


REZ- DT2 


System of equations written in matrix notation 


-1 

X := M -R X contains the constants of integration 

X={Clf ,Clc,C2c,Clm / C2m, ez } 
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set-up increments in radial location 


N := 10 i := 0 ..N rc := a rc := b rm := b 

0 10 0 


fiber compliant layer matrix 

a b - a c - b 

rf := — i rc := i + a rm := i + b 

i N i N i N 


FIBER 

radial displacment 

1 + VF DT 

Uf := AF rf + X • rf 

i 1 - VF 2.0 i 0 i 

radial stress 

DT 

SRf := -A + B X + B VF X 

2.0 0 5 


hoop stress 
-A 

STf := DT + B X + B VF X eqn (A— 3 ) 

2.0 0 5 

longitudinal stress 


SZf := -A- DT + B-2.0 VF X + B- (1 - VF) • X eqn (A-4) 

0 5 


eqn (A-l) 


eqn (A-2) 
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COMPLIANT LAYER 


radial displacment 


1 + VC DT 

Uc := AC 

i 1 - VC 2.0 


rc - 


a 

rc 


+ X • rc + — 
1 i rc 


radial stress 


DT 

SRc := -CC- 

i 2.0 


■ 


2' 


X 

T A 

a 


+ CD- 

1 

1*0 "" 

rc 


1 - 2 - VC 



i_ 










m 


+ CDM- VC- X 


hoop stress 



" 


2' 


DT 


a 



STc := CC 

1.0 - 



+ CD- 

i 2.0 


rc 




- 


- 



X X 

1 2 
+ 


1 - 2- VC 


[ rc J . 


- CC DT + CDM- VC- X 


longitudinal stress 


SZc := -CC- DT + CD- 2.0 


VC 

X + CDM- (1 - VC) • X 

(1 - 2 VC) 1 5 
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MATRIX 


radial displacement 


1 + VM DT 

Um := AM 

i 1 - VM 2.0 


rm - — 
i rm 


+ X - rm + — 
3 i rm 


radial stress 


SRm := -c- 
i 





2 


X 

X 

DT 

l.o - 

' b ' 


+ D- 

3 

4 






o 

• 

CM 


rm 



1 - 2-VM 

2 







Trm ] 







L iJ J 


+ DM VM- X 


hoop stress 





2 


DT 


' b ' 



STm := C- 

1.0 - 



+ D- 

i 2.0 


rm 




- 


- 



X X 

3 4 

+ 


1 - 2 - VM 


["*]' 


- C- DT + DM- VM- X 


longitudinal stress 

VM 

SZm := -C-DT + D-2.0 X + DM (1 - VM) - X 

(1 - 2-VM) 3 5 
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APPENDIX B 


ELASTIC 

Assuming temperature— dependent material properties (i.e., E(T) and a(T)) the 
stress— strain relations in equation (1) become: 

Aa r = _ E ^ { (l-i/)Ae r + v(Ae & + Ae z ) - (1 +i/)oAT} 

+ _ (1— ^)e + v(e $ +e) - (l+i/)o(T-T )} 

( 1-2 v){l + v) 

- E Aa(T— T ) (Bl) 

C 1—2*/) 


Aa fl = { (l-i/lAe* + v{Ae +AeJ - (l+i/)aAT} 

0 ( i-2 o ( irry 9 r z 

+ AE { (l-^)e„ + K»,+ e z ) - (l+^oCr-T )} 

( 1-2 1 / ) (1 + 1 /) 

- E Aa(T-T ) (B2) 

( 1 - 2-1 


Acr z = v{A a r + Aa^) + E(Ae z - aAT) + AE(e z - a(T— T Q )) 

-A a E(T— T q ) (B3) 

where A indicates an incremental change in the variable for a given increment in 
temperature (AT). Upon examining the above equations, one can sense the importance of 
including the effect of temperature— dependent material properties when AE and A a are 
large; that is, when E and a are strongly temperature— dependent. 

ELASTIC - PLASTIC 

Similarly, if we assume temperature— dependent elastic— plastic behavior, the 
multi— axial incremental stress— strain equations can be written in index notation (ref. 16) 
as follows: 


ij 


= V T > 


•kl 


+ C ijT 


(B4) 


where 


3D 


V 


i jkl 


(T) 


dT 


E kl 
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Decomposing the increment of total strain into an additive sum of elastic, plastic, and 
thermal contributions 


e ..= e e .+ ;?.+ e th 

ij 'i U 


'e ' p th 


(B5) 


and assuming a Von Mises yield condition, the increment in plastic strain is defined as 


( ij= dAS ij 


where S , j D. j k , ( e k ,-i ‘ f ) + (S ij C jj -fYft)T 


2 v dY 


dA = 


4 Y dY 
^ 5E P 


+ ^mn^mnpq ^pq 


and taking 


e 11 } = a • -T 
U iJ 


with 


Ojj = «fY + 1 J (T-T ) 

IJ 1 J firp ° 


da ■ ■ 


dT 


Given the above, we obtain the following expression for the stress rate, 


(B6) 


where 


a-. = L..,,(T) e,, +H..T 
ij ljkD > kl ij 

D- • S S D , , 
T _ n l jmn mn pq pqkl 

ijkl - U ijkl " t 


H.. = C..-T,.. u „ u - D *jkl S kl( S pq C pq - f Y IT ) 

ij ij ijkl kl g 


and 


^ ^ Y "cTe p + ^ mn ^ >mn P c l ^PQ 


where Y and E p are the equivalent yield stress and plastic strain, respectively. Again it is 


32 


apparent that if and Y are strongly dependent upon temperature, their 

influence should be taken into account. 


TABLE I . - BOUNDARY CONDITIONS FOR A 
CONCENTRIC CYLINDER MODEL WITH 
THREE CONSTITUENTS, FREE SURFACE 
CONSTRAINT, AND GENERALIZED 
PLANE STRAIN END CONSTRAINTS 


Locat ion , 
r 

Condi t ion 

0 

f 

u (0) = f ini te value 

a 

u f (a) = u c (a) 

a 

a p( a ) = <^(a) 

b 

u C (b) =■ u m (b) 

b 

°>) - aj(b) 

c 

o m T (c) = 0 


and equilibrium in the z 
direction , i ,e . , 


f r b c 

o tit dr + o_ irr dr 
J n Z J Z 


a m irr dr » 0 
b z 

where superscripts f, c, and m 
represent the displacement and 
stresses corresponding to the 
fiber, compliant layer, and matrix 
cylinders , respectively. 
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TABLE II. - COMPARISON OF SURFACE AND END BOUNDARY CONDITIONS 
SHOWN IN FIGURE 3, CONSIDERING TEMPERATURE- I NDEPENDENT 
ANALYTICAL SOLUTION PROCEDURE 


End 

condi lions 

Location, 

r 


Surface constraint 



Free boundary 

Homogenized 





Stress, ksi 





a r 

a e 

°z 

°r 

<*e 

a z 

Plane strain 

Fiber 

-32.8 

-32.8 

145.6 




Uz - O) 

Matrix: 

r - a 

-32.8 

76,5 

125.4 





r - b 

-23.4 

67.0 

125.4 





r * c 

0 

45,1 

125.4 





Fiber 

-29.5 

-29.5 

147.2 





Compliant layer: 
r - a 

-29.5 

111.0 

135.2 





r - b 

-17.3 

98.8 

135.2 





Matrix: 

r - b 

-17.3 

49.7 

122.4 





r = c 

. 

0 

33.4 

122.4 




Generalized 

_ — , — - 

Fiber 

-32.6 

-32.6 

-90.4 

-37.5 

-37.5 

-125.4 

plane strain 
(e z * 0) 

Matrix: 

r - a 

-32.6 

76.0 

60.1 

-37.5 

73 .2 

49,2 

r - b 

-23,2 

66.5 

60.1 

-27.9 

63.7 

49.2 


r * c 

0 

44.7 

60.1 

-4.3 

40.0 

49.2 


Fiber 

-29.3 

-29.3 

-93.6 

-34.9 

-34.9 

-129.6 


Compliant layer: 
r * a 

-29.3 

110.7 

101.9 

-34,9 

108.2 

95.2 


r = b 

-17.1 

98.6 

101.9 

-22.6 

95.9 

95.2 


Matrix: 

r - b 

-17.1 

49.2 

55.8 

-22.6 

45.6 

44.4 


r - c 

0 

33.1 

55.8 

-5.0 

28.0 

44.4 
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TABLE III. - GEOMETRY, LOAD, AND MATERIAL PARAMETERS FOR 

ANALYSES ASSOCIATED WITH TABLE II 

[Radii: a = 0.632; b = 0.6952; c = 1.0. 

Compliant layer thickness, t/a =0.1; 
temperature difference, AT = -1425 °F ; matrix 

yield at room temperature, = 53.89 k s i . ] 


Material property 

Fiber 

Compl iant 
layer 

Matrix 

Coefficient of 

thermal expansion, a, 
in . /in . /°F 

1 .96xl0~ 6 

lO.OxlO" 6 

S.OxlO" 6 

Modulus of elasticity, 
E, Msi 

58 

8.0 

16 

Poisson ratio, v 

0.25 

0.26 

0.26 


TABLE IV. - COMPARISON BETWEEN ASSUMPTION OF COUPLING AND DECOUPLING LONGITUDINAL STRESS TO [N-PLANE STRESSES 
FOR THREE COMPLIANT LAYER THICKNESSES (i.o., I/a = 0, 0.12, AND 0.2.) subjected to AT - 1472 °F 


Matrix 
( inner 
radius) 

No compl iant 1 ayer 
a * b - 1 .0 ; c 

(t/a = 0): 
1 .58 

Compliant layer (t/a - 0.12): 
a * 1.0; b - 1.12; c = 1.77 

Comp 1 i 
a - 1 .( 

int layer (t/ 
3 ; b - 1.2; 

'a * 0.2): 
* 1 .897 

Decoupled 

Coup 1 ed 

Oi f ference , a 
percent 

Decoupled 

Coupled 

Di [Terence , a 
percent 

Decoupled 

Coupled 

Di f Terence , a 
pnrcen t 

°r 

-20,7 

-29.77 

3.6 

-22.77 

-26 .85 

15 

-20.0 

-25.32 

21 


66.8 

69.57 

4 

53.2 

62.7 

16 

46.7 

59.12 

21 

°z 

46.9 

59.53 

21 .2 

41 .02 

55.58 

25 

38.1 

53.2 

28 

e z 

-0,005167 

-0.005204 

0.7 

-0.005423 

-0.005457 

0.6 

-0.00563 

-0.005617 

0.2 


difference is found by ( O couplod a decoup]ed)/ CT couplt'd 


Material property 

Fiber 

Compl i ant 
layer 

Matrix 

Coefficient of 

thermal expansion, x, 
in./in./°F 

2 .72xl0“ 6 

4 .38x 10~ 6 

6.5xt0" 6 

Modulus of elasticity, 
E, Msi 

62 

10 

10.9 

Poisson ratio, v 

0.25 

0.30 

0.30 
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TABLE V. - EXPERIMENTALLY OBTAINED (SEE REFS. 11 AND 12) TEMPERATURE-DEPENDENT 

FIBER AND MATRIX MATERIAL PROPERTIES 
(a) Fiber (SiC; SCS-6) 


Material property 

Temperature, °F 


77 

214 

398 

571 

753 

Coefficient of 

thermal expansion, a, 
in./in./°F 

1. 96xl0- 6 

1.98x10-6 

2.01x10-6 

2.07x10-6 

2.15x10-6 

Modulus of elasticity, 
E, Msi 

58 

58 

58 

58 

58 

Poisson ratio, v 

0.25 

0.25 

0.25 

0.25 

0.25 


Material property 


Temperature, ° 

F 



932 

1109 

1296 

1472 

1652 

Coefficient of 

thermal expansion, a, 
in,/in./°F 

2 ,24xl0 -6 

2.33x10-6 

2.42xl0- 6 

2.48xl0‘ 6 

2.55xl0" 6 

Modulus of elasticity, 
E , Msi 

58 

58 

58 

58 

58 

Poisson ratio, v 

0.25 

0.25 

0.25 

0.25 

0.25 


(b) Matrix (Ti— 24AI — 1 iNb ) 


Material property 



Temperature, °F 




75 


797 

1112 

1202 

1500 

Coefficient of 

thermal expansion, a, 
in./in,/°F 

5.0x10-6 

5.2x10-6 

5.7x10-6 

5.85x10-6 

5.9x10-6 

6.15xl0- 6 

Modulus of elasticity, 
E, Msi 

16.0 

14.5 

11.0 

12.5 

9.89 

6.2 

Stress, ay, ksi 

53.89 

59.0 

53.7 

42.2 

39.1 

24.0 

Hardening slope, 
H, Msi 

3.333 

0.441 

0.322 

0.187 

0.097 

0 

Poisson ratio, v 

0.26 

0.26 

0.26 

0.26 

0.26 

0.26 
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TABLE VI. - COMPARISON BETWEEN STRESS RESULTS USING TEMPERATURE-DEPENDENT 
AND TEMPERATURE- INDEPENDENT MATERIAL PROPERTIES 


[t/a = O.l; E C /E m = Oy/Oy = H C /H m = 0.5. ] 


Location 

Tempera ture-dependen t 
proper t i es 

Temperature- independent 
proper t ies 

Di f ference, 3 
percent 

a r 

a e 

°z 

a r 


a z 

r 

0 

z 

Fiber 

-33.0 

-33.0 

-105.0 

-29.4 

-29.4 

-93.7 

10.9 

10.9 

10 

.7 

Compliant layer: 











Inner radius 

-30.6 

122.1 

114 .2 

-27.3 

109.1 

102.0 

10.8 

10.6 



(r = 0.636) 











Outer radius 

-20.4 

112 .0 

114.2 

-18.2 

100.0 

102.0 

10.7 

10.7 



(r = 0.683) 











Matrix: 











Inner radius 

-16.0 

52.2 

114.2 

-14.3 

46.6 

55.9 

10.6 

10.7 



(r = 0.721) 











Outer radius 

-1 .0 

37.2 

62 .6 

-.9 

33.2 

55.9 

10 

10.7 



(r = 0.964) 












aDifference found by [()dependent _ ( ) independent } ]/( jdependent . 


TABLE VII. - COMPARISON OF ELASTIC AND F.I.ASTIC-PLASTIC STRESS ANALYSIS RESULTS 
[t/a - 0.1; F. c /E m = - H c /ll nl - 0.5. 1 
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TABLE VIII. - CASES ("UNKNOWN” 
COMPLIANT LAYER MATERIAL 
PROPERTIES) INVESTIGATED 
RELATIVE TO MATRIX 
MATERIAL PROPERTIES 



TABLE IX. - RADIAL RESIDUAL STRESSES AT INDICATED RADIAL LOCATIONS 
FOR VARIOUS CASES DESCRIBED IN TABLE VIII 

[t/a = 0.1 and a C /a m = 2.O.] 


Case 

Locat ion 


Fiber 

Compliant layer 

Matrix 



Inner radius 

Outer radius 

Inner radius 

Outer 

radius 




Radial stress, 

a r , ksi 



1 

-19.0 

-18.1 

-14,4 

-11.7 

-0 

.5 

2 

-24.4 

-22.7 

-15.0 

-11,9 

- 

.6 

3 

-20.2 

-19,1 

-14,2 

-11.4 

- 

.5 

4 

-17.5 

-16.8 

-13.8 

-11.3 

- 

.5 

5 

-23.3 

-21.6 

-14.8 

-11.7 

- 

.6 

6 

-21.9 

-20.5 

-14.7 

-11.7 

- 

.55 

7 

-19.1 

-18.2 

-14.0 

-11.4 

- 

.5 

8 

-24.2 

-22.5 

-15.2 

-12.0 

- 

.6 

9 

-20.0 

-18.9 

-14.5 

-11.7 

- 

.5 

10 

-19.0 

-18.1 

-14.2 

-11.5 



11 

-18.4 

-17.6 

-14.1 

-11.4 



12 

-19.1 

-18.2 

-14.2 

-11.5 



13 

-20.4 

-19.3 

-14.6 i 

-11.8 



14 

-21.7 

-20.5 

-15.0 

-12.0 

-.6 

15 

-21.3 

-20.1 

-14.9 

-11.9 

- 

5 

16 

-15.6 

-15.2 

-13.3 

-10.9 

- 

5 

17 

-19.1 

-18.1 

-13.9 

-11.3 

— 

5 
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TABLE X. - CIRCUMFERENTIAL RESIDUAL STRESSES AT INDICATED RADIAL 
LOCATIONS FOR VARIOUS CASES DESCRIBED IN TABLE VIII 
[l/a - 0.1 and a c /a ni « 2.0. ] 


Case 

Locat ion 


Fiber 

Compliant layer 

Matrix 



Inner radius 

Outer radius 

Inner radius 

Outer radius 


Circumferential stress, cq, ksi 

1 

-19.0 

36.8 

35.3 

35.2 

28.5 

2 

-24.4 

89.7 

86.8 

35.5 

28.9 

3 

-20.2 

52 .1 

52.5 

34.4 

27.7 

4 

-17.5 

27.0 

26.5 

34.1 

27.3 

5 

-23.3 

79.7 

76.4 

35.2 

28.6 

6 

-21.9 

66.3 

62.7 

35.2 

28.6 

7 

-19.1 

42 .4 

40.8 

34.4 

27.6 

8 

-24.2 

85.2 

82.3 

35.9 

29.4 

9 

-20.0 

46.1 

44.8 

35.2 

28.5 

10 

-19.0 

38.0 

38.0 

34.7 

28.1 

11 

-18.4 

34.5 

32.7 

34.6 

27.8 

12 

-19.1 

39.5 

37.9 

34.8 

28.1 

13 

-20.4 

49.5 

48.3 

35.3 

28.8 

14 

-21.7 

59.5 

58.7 

35.8 

29.4 

15 

-21.3 

57,3 

54.9 

35.8 

29.2 

16 

-15.6 

11.9 

12.8 

33.2 

26.4 

17 

-19.1 

43.7 

42.1 

34.1 

27.3 


TABLE XI. - LONGITUDINAL RESIDUAL STRESSES AT INDICATED RADIAL LOCATIONS 
FOR VARIOUS CASES DESCRIBED IN TABLE VIII 
[ t/a = 0.1 and o^/a 01 = 2.0. ] 
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TABLE XII. - ORDERING OF CIRCUMFERENTIAL AND LONGITUDINAL COMPLIANT LAYER STRESS 
COMPONENTS FROM MINIMUM TO MAXIMUM VALUES 
[t/a = 0.1 and cx c /a m = 2.0. ] 


Case 

Inner 

Outer 

Inner 

Out 

er 

Cor responding materi 

al properties 


rad 

ius 

radius 

radius 

rad 

i us 








Ci rcumfi 

arent ial 

Longi 

tudinal 

H c /H m 

c / m 

V a y 

c, m 

V°y 

H C 

/H m 

E C /E m 



stress , 

stress , 










ksi 

a z » 
ksi 


( 

i) 

(2) 

(3) 

(4) 

16 

11 

.9 

12.8 

11.5 

14 

.3 

0 

.5 

0.5 

0 


0.666 

4 

27 

.0 

26.5 

25.8 

28 

.3 


.75 

.5 


.25 

.666 

11 

34 

.5 

32.7 

32.6 

34 

.3 


.875 

.375 


.5 

.666 

1 

36 

.8 

35.3 

34.9 

37 

.1 

1 

.0 

.5 


. 5 

.5 

10 

38 

,0 

38. 0 a 

36.5 

40 

.3 a 



.75 


.25 

.666 

12 

39 

.5 

37. 9 a 

37.4 

39 

.8 a 



.5 


.5 

.666 

7 

42 

.4 

40.8 

40.3 

42 

.9 



.5 



1.0 

17 

43 

7 

42.1 

41.7 

44 

.3 



.5 



1 .25 

9 

46 

1 

44.8 

43.9 

47 

.1 

1, 

,16 

.666 



.666 

13 

49 

5 

48.3 

47.2 

50 

7 

1 , 

,25 

.75 



.666 

3 

52, 

1 

52.5 

50.1 

55 

,5 

1 , 

25 

1.0 


25 

1.25 

15 

57, 

3 

54.9 

53.8 

56 

.8 

1 , 

,5 

.5 

1 

0 

.666 

14 

59, 

5 

58.7 

56.9 

61 

,6 

1 . 

,5 

1.0 


5 

.666 

6 

66. 

.3 

62.7 

62.2 

64, 

,5 

1 . 

5 

.5 

1 

0 

1 .0 

5 

79, 

,7 

76.4 

75.2 

78, 

,8 

1, 

75 

.75 



1.25 

8 

85. 

.2 

82.3 

80.8 

85, 

,3 

2. 

0 

1.0 



1.0 

2 

89. 

,7 

86.8 

84.9 

89. 

,7 

2. 

0 

1.0 



1.25 


a The only exceptions to the ordering. 
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TABLE XIII. - PURE-ELEMENT CANDIDATE COMPLIANT LAYER MATERIALS FOR 
SCS-6/Ti 3 AI + Nb COMPOSITE SYSTEM 


Element 

Coefficient of 
thermal expansion, 

ct/a" 1 

Yield stress, 

/ m 
a /o 

y y 

Modulus of, 
elast ici ty , 

E/E m 

Melting 

point, 

°F 

Comment 

Melting point above 1500 °F 

Europium 

2.88 


n ire; 

1519 



U . LOO 


Ytterbium 

2.78 

0.02 

.22 

1515 


Calcium* 

2,48 

.04 

.22 

1540 

Handleabi 1 i ty in air 

Manganese 

2.44 

.65 

1.44 

2273 


Si lver 

2.18 

.15 

.687 

1761 


Copper 

1 .84 

0.09 

1 .0 

1981 


Gold 



.725 

1965 


1 . 58 



Coba 1 1 * 

1.54 

2.04 

1.87 

2723 

Too high a yield 

Nickel 

1.48 

.16 

1.875 

2647 


Thorium* 

1.38 

.385 

.652 

3182 

Radioactive 

Iron 

1.3 

.16 

1.812 

2798 






2826 


Pal ladium 

1 .3 




Beryl 1 ium* 

1.28 

.5 

2.5 


Carcinogenic 

Melting point between 1000 and 1500 °F 

Plutonium* 

6.11 

0.74 

0.875 

1184 

Radioactive 

Magnesium 

3.01 

.24 

.375 

1202 


Aluminum 

2.62 

.05 

.56 

1220 



‘Unlikely candidate; see comment. 


TABLE XIV. - ALLOY CANDIDATE COMPLIANT LAYER MATERIALS 


FOR SCS-6/Ti 3 Al + Nb COMPOSITE SYSTEM 


A1 loy 

Coefficient of 
thermal expansion, 

<x/a m 

Yield stress, 

/ m 
a /a 

y y 

Modulus of, 
elastici ty , 

E/E m 

Melting 
point , 
°F 

Comment 

Yel low brass 

2 

.1 

0.33 

0.875 

1710 


Ni-Ag cast 

2 

.0 

.46 


1980 


(20 percent ) 







Red brass 

1 

.96 

.3 

.937 

1875 


Si I icon bronze 

1.9 

.41 

.937 

1865 

— 

Phosphor bronze 

1 

.88 

.37 

.937 

1920 


Stainless steel 

1 

.86 

.55 


2550 


321, 347 







A1 bronze 

1 .84 

.46 

.937 

1900 


Incoloy 800 

1 

58 

.67 

1.77 

2475 


Incoloy 801* 

1 

56 

1.04 

1 .88 

2475 

Too high a yield 

Incoloy 802* 



.78 

1,85 

2450 

Too high a yield 

Durimet T 



.65 

1.42 

— 


Monel 



.65 

1.625 

2370 


N-155* 

1.55 

1.08 

1.83 

2325 

Too high a yield 

Inconel 690* 

1 

5 

.83 

1 .9 

2450 

Too high a yield 

Hastel loy-X* 

1.54 

.96 

1.78 

2300 

Too high a yield 


‘Unlikely, candidate; see comment. 
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ORIGINAL' PAGE 

BLACK AND WHITE PHOTOGRAPH 



(a) As fabricated. (b) After 1000 thermal cycles. 

Figure 1.— Experimentally observed radial cracking in SICAV24AI-1 INb after fabrication and after an additional 1000 thermal cycles, 
taken from reference 1 . 



(a) Possible internal crack orientations. (b) Unit cell and coordinate system. 

Figure 2. — Definition of concentric cylinder model and possible internal crack orientations. 
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l-Fp'Vj 



End constraints 



mmzm 



y/zy//////////A 

E Z =0 




Free Free Free 

(a) Generalized plane strain. (b) Plane strain. (c) Homogenized. 

Figure 3. — Various assumed end and surface boundary conditions. 




Figure 4. — Finite element concentric cylinder idealization with 
a compliant layer thickness t/a of 0.05. 


Figure 5.— Finite element concentric cylinder idealization with 
a compliant layer thickness t/a of 0.10. 
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X Arithmetic mean location 
• Integration point 



2 


Figure 6. -—Illustration of MARC 8-node brick element and 
8-point Gauss quadrature integration scheme. 



r, in. 

Figure 8.— Comparison between finite element and ana- 
lytical circumferential stress versus radial location. 



r, in. 


Figure 7. — Comparison between finite element and ana- 
lytical radial stress versus radal location. 



r, in. 

Figure 9. — Comparison between finite element and ana- 
lytical longitudinal stress versus radial location 
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(a) Fiber/matrix. (c) Proposed solution. 


Figure 10.— Description of thermal mismatch problem and compliant/ 
compensating layer solution. 



Figure 1 2 — Effect of varying CTE ratio a c /a m 
and thickness t/a of compliant layer on elastic 
longitudinal stress distribution. 



Figure 1 1.— Effect of varying CTE ratio a c / a m and 
thickness t/a of compliant layer on elastic circum- 
ferential stress distribution. 



Figure 13.— Effect of varying stiff- 
ness of compliant layer (t/a =*0.1) 
on elastic circumferential stress 
distribution in fiber and matrix. 
Inner radius. 
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Figure 18 —Effect of varying CTE ratio a 9 /a m and thickness 
t/a of compliant layer on elastic-plastic radial stress distribution 
in compliant layer. 





Figure 19. — Effect of varying compliant layer hardening slope 
H c /H m on circumferential and longitudinal (elastic-plastic) 
stress distribution in compliant layer. 



Figure 20.— Effect of varying compliant layer yield point Oy/a™ 
on circumferential and longitudinal (elastic-plastic) stress distri- 
bution in compliant layer. 
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Figure 21 . — Effect of varying compliant layer stiffness E° /E m 
on circumferential and longitudinal (elastic-plastic) stress distri- 
bution in compliant layer. 


IR 



Figure 22. — Longitudinal matrix stress distribution, sum- 
marized using inner and outer radial location values, 
versus normalized CTE ratio a c /a m for different 
compliant layer thicknesses (i.e., t/a = 0,05, 0.1, and 
0.2). 


48 


Report Documentation Page 

National Aeronautics and r ^ 

Space Administration 

t Report No. 2. Government Accession No. 

NASA TM- 103204 
4. Title and Subtitle 

Elastic/Plastic Analyses of Advanced Composites Investigating the Use of 
the Compliant Layer Concept in Reducing Residual Stresses Resulting 
From Processing 

7. Author(s) 

Steven M. Arnold, Vinod K. Arya, and Matthew E. Melis 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


Recipient’s Catalog No. 


Report Date 
September 1990 
Performing Organization Code 


Performing Organization Report No. 

E-5661 

Work Unit No. 

510-01-01 
Contract or Grant No. 


Type of Report and Period Covered 
Technical Memorandum 

Sponsoring Agency Code 


15. Supplementary Notes 

Steven M. Arnold, Lewis Research Center; Vinod K. Arya, University of Toledo and National Research Council 
NASA Research Associate; Matthew E. Melis, Lewis Research Center. 


16. Abstract 

High residual stresses within metal and intermetallic matrix composite systems can develop upon cooling from the 
processing temperature to room temperature due to the coefficient of thermal expansion (CTE) mismatch between 
the fiber and matrix. As a result, within certain composite systems, radial, circumferential, and/or longitudinal 
cracks have been observed to form at the fiber-matrix interface region. The compliant layer concept (insertion of 
a compensating interface material between the fiber and matrix) has been proposed to reduce or eliminate the 
residual stress buildup during cooling and thus minimize cracking. The present study investigates both elastically 
and elastic-plastically the viability of the proposed compliant layer concept. A detailed parametric study was 
conducted utilizing a unit cell model consisting of three concentric cylinders to determine the required character 
(i.e., thickness and mechanical properties) of the compliant layer as well as its applicability. The unknown 
compliant layer mechanical properties were expressed as ratios of the corresponding temperature dependent 
Ti-24Al-llNb (a/o) matrix properties. The fiber properties taken were those corresponding to SCS-6 (SiC). 
Results indicate that the compliant layer can be used to reduce, if not eliminate, radial and circumferential 
residual stresses within the fiber and matrix and therefore also reduce or eliminate the radial cracking. However, 
with this decrease in in-plane stresses, one obtains an increase in longitudinal stress, thus potentially initiating 
longitudinal cracking. Guidelines are given for the selection of a specific compliant material, given a perfectly 
bonded system. 


17. Key Words (Suggested by Author(s)) 

Residual stress; Thermal mismatch; Elastic; Plastic; 
Compliant layer; Finite element 


18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 24 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No. of pages 

50 


22. Price* 

A03 


NASA FORM 1626 OCT £ 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 




